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A general formulation of translationally invariant, parametrically correlated random matrix ensem- 
bles, is used to classify universality in correlation functions. Surprisingly, the range of possible 
physical systems is bounded, and can be labeled by a parameter a £ (0, 2], in a manner analogous 
to Levy diffusion. Universality is obtained after scaling by the (anomalous) diffusion constant D a 
(the usual scaling is divergent for a < 2). For each a, correlation functions are universal, and 
distinct. The previous results in the literature correspond to the limiting case of superdiffusion, 
a = 2. 
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It is becoming increasingly possible to explore the 
^^physical properties of complex quantum systems with the 
'tools offered by Random Matrix Theory (RMT). One 
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^ of the main lessons learned from systematic studies of 
Oc omplex quantum many-body systems and simple quan- 
.^^tum systems with a classical chaotic limit, is that they 
display the same kind of spectral universal fluctuations, 
^^depending only on the type of symmetry that they ex- 
hibit, e.g. orthogonal, unitary or symplectic [1]. Until 
very recently, however, very little was know about chaotic 
^systems which depend on an external parameter x, and 
f~^how the properties of such systems are correlated in that 
^fparameter. While the parameter is often viewed as an 
^^external quantity, such as an electric or magnetic field, it 
also important in the study many-body systems with 
^C^slow and fast modes, such as atomic nuclei. In that case, 
Git corresponds to the slow variables in the adiabatic limit 
the time-evolution of the many-body system. One 
I of the remarkable properties which has emerged is that, 
Ounder appropriate scaling, parametric-dependent correla- 
ion functions are universal as well. One of the first works 
^addressing such question was a study of level curvature 
distributions [2]. A significant step towards the under- 
standing of parametric correlations was recently achieved 
by Szafer, Simons and Altshuler [3,4] in a sequence of pa- 
pers. They were able to compute analytically the density- 
density and level velocity correlators in a particular ran- 
dom matrix model realization, which is very successful 
in describing different chaotic and disordered systems. 
Correlation functions were found to fall upon universal 
curves, seemingly independent of the underlying proper- 
ties of the studied systems. In this letter, we show that 
there are in fact a continuous number of universal curves 
for each type of observable, and that those known up to 
now are a just a limiting case. Furthermore, the quantity 
which is used to rescale the parameter to obtain univer- 
sality, is generally divergent, requiring the introduction 
of a more general scaling. 

The framework which we develop here is based on 
assuming that we have an ensemble of random matrix 



Hamiltonians H(x), which depends on an external pa- 
rameter x. At each value of x, the matrix has the same 
spectral properties, pertaining to the Gaussian Orthogo- 
nal Ensemble (GOE) (or equivalently the Gaussian Uni- 
tary (GUE) or Symplectic (GSE) Ensembles; we will con- 
sider here only the GOE case, as the others can be ob- 
tained in a straightforward manner). Furthermore, we 
ask that the system have translational invariance, in the 
sense that correlation functions at two different values of 
the external parameter x, y depend only on their separa- 
tion \x — y\. Hence, our starting point is the assumption 
that the covariance of H is translationally invariant: 
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H ij (x)H kl (y) = -F(x- 



(l) 



Here F(x) = F(—x) defines the rate of decorrelation of 
the system as one varies the parameter. There are two 
independent formulations we use to realize such a H — H 
correlation. 

(I)Stochastic Integral: In this first approach, we consider 
the explicit construction of H(x) as the most general lin- 
ear combination of independent, uncorrelated GOE ma- 
trices: 



H ij (x) = J dyf(x-y)V tJ (y). (2) 
Here Vij(y) is uncorrelated white noise at every i,j,y: 



Vij(y) = 0, Vij(y)Vu(y') = 6(y - y')(h ik bji + 6 jk 6u) . 

(3) 

The stochastic integral (2) defines a theory such as (1) 
through the convolution 



F(x -y)= I dzf(x - z)f(y - z) . 



(4) 



As H(x) must be a member of the ensemble for every x, 
we see that f(x) must be a real valued function. 
(II) Measure: Alternately, we can derive a generic covari- 
ance (1) from the measure: 
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P{H)VH oc Y[ dHij(x) x (5) 

xexp (-^ j j dxdyTr[H(x)K(x-y)H(y)]j . 

Since the measure is quadratic in H , the correlation func- 
tion F(x) is just the functional inverse of the operator K: 

j dz K(x- Z )F( Z -y) = 6(x-y). (6) 

The physical relevance of random matrix modelling 
is directly related with the identification of the scaling 
parameter with physical quantities. At the RMT level 
though, universality of correlation functions is typically 
obtained in the large N limit by scaling the parameter 
x with some power of N , where N is the dimension of 
H . As a consequence, the understanding of the short dis- 
tance behavior of F(x — y) plays a central role. To explore 
this, we choose a natural, reasonable and fairly generic 
parametrization of the entire range of possible leading 
order short distance behaviors of F(x — y) by choosing 
the function of the form: 

F(x-y) = e-\ x -»\ a ~l-\x-y\ a + ... . (7) 

Here a is an arbitrary real number. Surprisingly, we will 
see that the range of a cannot be arbitrary, but must be 
bounded. (As a result, parametric correlations are not 
generally related by a rescaling.) If we chose different 
realizations with the same small distance behavior, such 
as F(x) = 1/(1 + x a ), cos(x a l 2 ), and so forth, the con- 
clusions we find are identical. In this sense the specific 
choice of F is not important. 

Consider first the range of a for which the function 
K(x — y) defines a probability distribution. The measure, 
Eq. (5), is normalizable providing that for each pair i, j, 
the quantity 

j dx j dyH ij (x)H ji (y)K(x - y) > , (8) 

so that K must be a positive definite operator. Since 
F is the functional inverse of K, F must also be pos- 
itive definite. The properties of positive definite func- 
tions were systematically studied by Bochner [5]. Using 
the theorems of Bochner on positive definite functions 
[5], the only range of a in which F(x) is positive definite 
is < a < 2. Hence the measure is only defined for 
a£(0,2]. 

The conclusion is the same from method (I): F(x) can 
only be factorized into real functions / in Eq. (4) when 
F is a positive definite function. This is most easily seen 
in the special case when / itself is symmetric, then 

f(x ~ v) = h f dk y^) cos (M x - y)) > ( 9 ) 



where F(k) is the fourier transform of F(x — y). Again 
from Bochner we know that: the fourier transform F(k) 
is a positive function, if and only if F(x — y) is positive 
definite. Hence Eq. (9) is a well defined construction of 
f(x — y) for a £ (0, 2]. As a consequence, it is interesting 
to note that one cannot choose a general covariance (1) 
and find a construction (2): / can only be defined consis- 
tently when a £ (0, 2]. If F is not positive definite, then 
/ is in general complex, and H (x) is no longer a member 
of the chosen random matrix ensemble. 

We have explicitly constructed parametric Hamiltoni- 
ans H(x) for a £ (0, 2] by factorizing the positive definite 
function (7). In Fig. 1, we show a range of measured co- 
variances of H , using matrices of dimension N = 100, for 
the values of a = i,l,|,2 in Eq. (7). H(x) was con- 
structed over a range of x for each case, and F(x — y) 
was measured by averaging over x and y, the quantity: 

F ( x - y ^ = N( N-i) T, H 'i( x ) H 'i(y) ■ ( 10 ) 

The results in Fig. 1 are representative points, the grid 
size being finer, and are in agreement with the expected 
value of a. (For instance, a fit to the a = 3/2 curve 
yields an exponent of a = 1.49(1).) 
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FIG. 1. Computed covariance F(x — y) = H(x)H(y) from 
simulations at a = 1/2,1,3/2,2 in (7). The solid line is the 
exact result and the boxes are the results of the simulation. 
The x axis is rescaled in the figure in order to show the short 
range behavior. The points are shown are representative; the 
actual grid size is much smaller. 

To explore the influence of a on the form of univer- 
sal curves, we first consider the short distance behavior 
of the Hamiltonian as manifested in various observables. 
Using perturbation theory to explore the short distance 
properties, the energy and wavefunction overlaps, to sec- 
ond order, are 
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8E n (x)= E n (x') - E n {x) = (11) 
= W ""+L E _ E +- 

i^ n (,oi^w) i 2 =i-E (Mi + - (i2) 

Using the ensemble average defined by the correlation 
(1), the short distance behavior (7), and that we have 
translational invariance, it follows from Dyson [6] that 
(8E n (x)) 2 ~ 8x a , and 1 - |(*(a;')l*(*)}| 2 ~ 6x a . The 
case for a = 1 corresponds to Dyson's Brownian motion 
model for random matrices [6]. One immediate obser- 
vation is that not only do these quantities decorrelate 
differently for each a, but C(0), which is obtained from 

C(x - x 1 ) = C(6x) = 1F0*0WK) for 6x ~> °> d i ver g es 
as 1 /<5a? 2— a . By measuring C(0) as a finite difference on 
a grid defined by the discretization of Eq. (2), the di- 
vergence will be in terms of the minimum grid spacing. 
Thus C(0) is only finite for the special case a = 2, and 
one cannot obtain universality by scaling the parameter 
x by \J C(0). So, as <5a? — ► 0: 

C (5a; a - 2 a > 2 

C(0) ~ < constant a = 2 . (13) 

( (l/^a?) 2 -« oo a < 2 

(Our simulations have confirmed this behavior for a < 2.) 
It is clear that one must define a more general level veloc- 
ity, and a new scaling. Because wavefunction correlations 
do not suffer from this behavior, they serve as a better 
measure of the anomalous character a of the universal 
behaviour. Eq. (13) provides another argument on the 
finite range of a. If one had a realization of a > 2, 
one would have both energy level correlations from (1) 
and a vanishing rms value of the energy level slope since 
C(0) = 0, which is not possible. 

We would like to view the processes classified by a in 
terms of its similarities with Levy or anomalous diffusion 
[7]. Anomalous diffusion of a particle in a disordered or 
chaotic environment generally leads to diffusion with a 
mean square deviation of the type 

(Q 2 (t)) = DP , (14) 

with y zfz I. When j = 1, one recovers Brownian diffu- 
sion. In a similar manner, a natural extension of Dyson's 
Brownian motion model is then to interpret the behavior 
in terms of anomalous diffusion in the parameter (ficti- 
tious time): 

jE 2 ~=DJx a . (15) 

Here D a is defined as the anomalous diffusion constant. 
One would be tempted to object against this interpreta- 
tion arguing that there is a fundamental difference be- 
tween Levy flights (14) and the parametric result (15): 



The former is a long-time behavior, while the latter is re- 
lated to short-distance in parameter space. This contra- 
diction is only apparent, since the character of the con- 
sidered gaussian ensembles magnifies enormously small 
perturbations which results in effectively scaling the x- 
parameter by some power of N . As a consequence, the 
"short distance" leading order behavior become the dom- 
inant feature as N — ► oo. In this sense, the large N limit 
makes Eq. (15) analogous to Levy diffusion. For a = 1, 
we recover Brownian diffusion, for a > 1 one has su- 
perdiffusion, while a < 1 gives subdiffusion. Hence, the 
proper scaling of the coordinate must be related to the 
anomalous diffusion constant, which relates the average 
drift in energy to parametric distance: 

x=[D a ]^ a x, D " = J^- ( 16 ) 

Then C(0) can be identified as the diffusion constant 
for the superdiffusive case a = 2: D2 = C(0), so that 
x = y/C(0)x, recovering the universal scaling derived by 
Szafer, Simons and Altshuler [3]. 

In Fig. 2 we show the universal curves for the wave- 
function overlap obtained for the processes of Fig. 1. 
When plotted as a function of the x = [DaY^x, the 
decorrelations are universal, and distinct for each a. Re- 
call that scaling by %J C(0) does not result in universality 
since C(0) is divergent for a < 2. It is worth noting that 
the functional form of the wavefunction correlations are 
reproduced quite well by 

|(*(*0IW» 1 + c| ~_, r (17) 

which is Lorentzian for a = 2. (Numerical details, such as 
the nature of the constant c, will be discussed elsewhere). 
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FIG. 2. Comparison of universal wavefunction overlaps for 
the different anomalous processes of Fig.l. Universality is ob- 
tained by rescaling with the diffusion constant: x = [Da] 1 ' !. 
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We only point out here that the overlap | {^(a?') | ^(a?)) | 2 
provides a better measure of the value of a for a partic- 
ular physical system, since it is not singular like C(x). 
Hence, in this way, one can fit the functional form and 
measure a. 

There are several ways one can define a better mea- 
sure of the level velocities. One is to define a generalized 
diffusion constant (an extension of the definition of C(x) 
to include the stochastic measure for the anomalous pro- 
cess): 



D(x - y) 



6E(x) 6E{y) 
(6x) a f 2 (6x) a f 2 ' 



D(0) = D a 



(18) 



which for a = 2, reduces to D(x) = C(x). This function 
does not diverge on small scales. The form of D(x) / D(0) 
for a = 1/2, 1, 3/2, 2 is shown in Fig. 3. For a = 2, one 
recovers the well known curve, but is otherwise different 
for different a. Another is to introduce the fractional 
derivative of order a of the curve E(x), denoted A", as 



a finite difference quotient, 
construction [8] 



In the Griinwald-Liouville 
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A:E(x) = Km ^ £(-)* ( I ) E(x - k6x) . (19) 

k = 



Then the natural definition for the scaling parameter 
is: D(x) = 
5(0) ± D Q 



is: D(x) = A"/ 2 E(x)A™/ 2 E(y). Notice that in general 



shown that the range of possible short scale behaviours 
of the H — H correlator is bounded, and can be assigned 
a value a £ (0,2]. For every value of a, there exists 
universal correlation functions. However, one of the sur- 
prising results is that these functions are distinct for each 
a. They are not simply related by rescaling, due to the 
finite range of a. The conventional universal scaling by 
^/C(0) was seen to diverge for all but the superdiffu- 
sive a = 2. However, the introduction of scaling by the 
anomalous diffusion constant, was shown to be 

well defined for each a, and result in universality for each 
a. While we have constructed model Hamiltonians for 
different values of a, it would be nice to find realizations 
in nature. One way to explore this is to measure wave- 
function correlations, since these are non-singular, and 
one can extract the anomalous value a which governs 
the process. Finally, we note that when the parameter 
x = x(i) is dynamic, the solution to the time-dependent 
Schrodinger equation for different observables depends 
on the short time behavior of the quantity 1 — F(i) [9], 
and we suspect that a similar classification exists for time 
dependent quantities such as occupation numbers. 

This work was supported in part by the Department 
of Energy Grant DE-FG02-91ER40608. C. H. L. is sup- 
ported by the National Science Foundation. 
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FIG. 3. The level velocity correlator D(x)/D(0) versus x 
for a = 1/2, 1, 3/2, 2. The results for a = 2 are identical to 
those of Szafer, Simons and Altshuler Ref.[3]. (The a = 1/2 
case is numerically problematic since the slope of F(x) at 
x = diverges. For this reason the curve is more representa- 
tive, having been smoothed with a spline.) 

In conclusion, we have endeavored to formulate the 
most general class of correlated random matrix Hamilto- 
nians which are translationally invariant. From this we 
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